Hopping between localized Floquet states in periodically driven quantum dots 
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The dynamic localization in energy space - suppression of the absorption of energy from an 
external microwave field due to quantum interference - was analyzed recently for a closed quantum 
dot in the absence of electron-electron interactions. Here a weak interaction is shown to lead to a 
finite absorption and heating, which may be viewed as hopping between localized Floquet states. The 
heating rate grows together with the electronic temperature, eventually destroying the localization. 
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Introduction. — The dynamic localization (DL) - 
suppression of the absorption of energy by a quantum sys- 
tem under a periodic perturbation - was studied exten- 
sively forquantum chaotic models @,E3- ^ n a recen t pub- 
lication |j| DL was shown to be possible in a solid-state 
system - an ac driven chaotic quantum dot (QD). Exper- 
iments on such systems started in the last few years 
which makes the DL in a QD more than just a theoretical 
peculiarity. Calculations of Ref. 0] employed the sim- 
plest random-matrix model for the dot, which neglects 
electron-electron interaction. The latter, however, is al- 
ways present in a QD, and studying its effect is the sub- 
ject of the present paper. Many-body effects on the DL 
have been studied previously for quite a different sys- 
tem 0, so this issue is also of a fundamental interest. 

Let the single-electron mean level spacing 5 in the dot 
(always assumed to be closed) be the smallest energy 
scale in the problem. Then, if an external time-dependent 
periodic perturbation with the frequency u is applied, the 
total electronic energy E in the dot (counted from that of 
the ground state) grows linearly with time as described 
by the Fermi Golden Rule: E(t) = Tuj 2 t/{26) = W t. 
The probability of each single-electron transition per unit 
time, denoted by 2T, measures the strength of the pertur- 
bation [tJ. The criterion of validity of the Fermi Golden 
Rule is T 3> 6, and T <C oj is also assumed (h = 1). 

This picture (hereafter referred to as Ohmic absorp- 
tion) is valid provided that each act of photon absorp- 
tion by an electron is independent of the previous ones; 
however, for a discrete energy spectrum this turns out 
not to be the case. After many transitions the absorp- 
tion rate decreases due to accumulation of the quantum 
interference correction Q , so that after a time t* ~ T/S 2 
the absorption is completely suppressed. This effect was 
named the dynamic localization in energy space. It is 
the consequence of level discreteness: at 6 — » it takes 
longer time for the DL to develop, and for the contin- 
uous spectrum there is no DL. It should be also noted 
that the problem of DL in the energy space turns out to 
be analogous to the Anderson localization problem in a 
one-dimensional disordered system in the real space . 

Conduction in a disordered sample in the regime of the 
strong Anderson localization occurs via thermally acti- 



vated hopping between localized states due to inelastic 
processes (typically, phonons) the activated temper- 
ature dependence being one of the main signatures of the 
localization. In the QD phonons are assumed to be frozen 
out by the external cryostat, hence the dominant role is 
played by electron-electron collisions. The latter are gov- 
erned by the electronic temperature T in the dot, which 
is different from the cryostat temperature. In the DL 
regime the (effective) electronic temperature ~ Tlu/5 
and the localization length in the energy space are in fact 
the same thing, which invalidates the familiar concept of 
the activated hopping conduction for the DL problem. 

Electron-electron collisions in an equilibrium QD were 
studied by Sivan, Imry, and Aronov (SIA) [j| (see also 
Ref. 0). The single-electron relaxation rate was found 
to be 7(T) ~ T 2 S/E 2 h , where E Th is the Thouless en- 
ergy. In the following, the hierarchy of scales 8 <C T <C 
to <C 2* -C Etk is assumed, Etk being thus the largest 
energy scale in the problem. For the strong DL to take 
place, the quantum correction should have enough time 
to develop, which is governed by a dimensionless param- 
eter u ee 7(T*)t* ~ (T/5) 3 {u/E Th ) 2 < 1. 

Qualitative picture. — The theory to be presented be- 
low results in the following qualitative picture. As the 
collisions are rare (ji* Cl), the electrons spend most of 
the time in the states localized in energy space, having 
definite phase relationships. When two electrons collide, 
the phase memory is lost for them, and their wave pack- 
ets start spreading along the energy axis. They localize 
again after the time ~ t*, in the meantime spreading 
by ~ T*. Thus, although an electron-electron collision, 
strictly speaking, by itself does not change the total en- 
ergy of the system, the subsequent ac driven dynamics 
leads to a change of the total energy of ~ T* per collision. 
The sign of this change is, however, arbitrary, because a 
periodic perturbation can equally cause transitions up 
and down the spectrum. Only the presence of the filled 
Fermi sea below (i. e., an energy gradient of the electronic 
distribution function) makes absorption the preferred di- 
rection, which means that if the electronic temperature 
T ^> T», the energy absorbed per collision is ~ T 2 jT 
rather than T*. The effective number of electrons that 
can participate in a collision is ~ T/5 (due to the degen- 
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FIG. 1: A sketch of the dependence of the dimensionless ab- 
sorption rate W(t)/Wo on time in the case of strong dynamic 
localization (it = 7(T»)£» <C 1)- Initially the absorption is 
Ohmic, at t ~ t* it is suppressed by the DL to a small resid- 
ual value W ~ uWo (note the algebraic smallness in contrast 
to an exponential one for hopping conduction in disordered 
media). After a long time t*/u = l/j(T*) the exponential 
heating becomes noticeable, and after a time (i*/it) ln(l/w) 
the strong DL is destroyed and W ~ Wo- 



erate Fermi statistics). During the time interval ~ I/7 
each of these electrons participates in one collision, so the 
total number of collisions per unit time is ~ (T/6)^(T). 
This gives the energy absorption rate: 



Trr dE 7? T , , 
W = — ~ — — 7 (T) 
dt T 8 



T Tj 

E Th 



l(T*)E, (1) 



which is the main result of the paper. The total energy 
(E ~ T 2 /S) is growing exponentially, however, with the 
rate being small in the parameter u. This picture is valid 
up to T such that 7(T) ~ 1/t*, when the strong DL is 
destroyed. One can see that at this temperature Eq. Q 
gives the Ohmic absorption rate Wo ~ Yuj 2 /S ~ wT,. 
This happens for an arbitrarily weak interaction, in con- 
trast to the case of a kicked nonlinear quantum rotor Q . 

Eq. JD allows for another simple interpretation. Each 
single electron is hit by another one with the average 
time interval 1/7. After a collision the electron spends 
the time i» <C 1/7 spreading its wave packet, and then 
it has to wait for the next collision. Thus the absorption 
rate of the whole system is given by the simple weighted 
average: W ~ Wot*-y(T), which is exactly Eq. Q. 

Floquet states. — It is convenient to start the for- 
mal part of the discussion from some general state- 
ments concerning a system whose hamiltonian H(t) is 
periodic in time with a period 27r/w. According to 
the Floquet theorem the eigenfunctions of the time- 
dependent Schrodinger equation have the form \ip a (t)) — 
e -i\ a t |^ Q ,(t)) ) where \ X a(t)) is a l so periodic in t, and 
the range of quasienergies is restricted to |A Q | < uj/2. 
One can show that |Xa(*)) a t any fixed t form an 
orthonormal and complete system: 

<Xa(t)\Xa'(t))=6 aa >, X)|Xa(*))(Xa(t)| = i, (2) 

a 

where 1 is the unit operator. 



If some perturbation H'(t) is added to the original 
hamiltonian H(t), the solution of the Schrodinger equa- 
tion can be sought in the form \i(j(t)) = J2 a C a (t) \Xa(t)), 
with the coefficients satisfying 

iC a (t) = \ a C a {t) + Y, (Xa(t)\H'(t) \ X a'{t)) C a ,(t). (3) 



Now suppose that some level A Q is resonant with a contin- 
uum of other states X a i , and H' does not depend on time. 
One can proceed analogously to the textbook derivation 
of the Fermi Golden Rule picking up the singular term 
in first-order perturbative expression for |C Q '(i)| 2 , and 
obtain the transition probability per unit time: 

dw a ^ a > = 2TT\M a > a \ 2 S(X a i - X a ) da', (4) 



M a , a = \ ( Xa '(t)\H'\ Xa (t)) 



uj dt 
~2w 



(5) 



-7T/u 



As A a ' may originate from the part of the spectrum of Hq 
which was higher or lower than A Q by several oj, Eq. lf4"|l 
describes photon-assisted transitions. 

Model for the dot. — The unperturbed single-particle 
hamiltonian Hq = — V 2 /(2m) + v(r) is the usual one, 
giving rise to the random-matrix-theory (RMT) spec- 
trum with the mean level spacing S at energies smaller 
than Eth ; 12j. A deterministic form of the perturbation 
V(t) can lead to the time-dependent RMT [13j, which is 
equivalent to another stochastic model: 

r<5 



V(t) = V(r) <j>{t) , V(r)V(r>) = —L d 6(r-r') (6) 

7T 



with V(r) random and independent of v(r), V(r) = 0, 
and <f)(t) being a given periodic function of time. Periodic 
functions are expanded in the Fourier series: 
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S — — OO 



(7) 

The harmonics 4> s are assumed to decrease faster than 
1/s 3 / 2 to give a finite absorption rate Q. Floquet wave 
functions satisfy the equation [Tlj : 

(Ao - S c)xi s) (r) + £ V(r)^ s , X { / } (r) = A Q xi s) (r). 

(8) 

No averaging over the disorder has been performed yet. 

Eq. (jSJ has a simple and illustrative interpretation. It 
describes stationary wave functions for a one-dimensional 
chain of coupled granules of the same shape, but with 
the energy of each granule s shifted by su>. The cou- 
pling broadens each level of an unperturbed granule by Y. 
Due to the condition Y ^> 5 the motion along the chain 
is diffusive with the dimensionless diffusion coefficient 
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(r/2) J2 S 1 0s | 2 . Due to the condition u > T the spec- 
tral correlations within the energy strip of the width T be- 
tween different granules can be neglected, and the gran- 
ules can be considered independent. Wave functions in 
such a system have been shown to be localized with the 
localization length £ ~ T/S ^> 1 14J . Neglecting fluctua- 
tions of the spectrum, one can label the localized states 
by the position s (the number of the segment of the size £) 
and the quasienergy A, introducing the mean level spac- 
ing 5(i = <5/£ ~ 1/t* within one segment: 



E 



ds 



u>/2 



-lu/2 



dX 



(9) 



Kinetic equation. — As the Floquet states provide a 
suitable single-particle basis, one can introduce their oc- 
cupation probabilities f a which are related to the actual 
energy distribution function f e via 



/ e = <^/ Q (xW|x«)<5(e-A a - SW ) 



(10) 



In the absence of electron-electron collisions f a does not 
depend on time, otherwise one can write down the colli- 
sion integral based on the Fermi Golden Rule QJ: 



dt 



2tt 



E 



M, 



CK3Q4 — *QiCK2 



I S(X ai +A a2 A Q3 A a4 ) [(1 f ctl )(l fot^fotzfotA faifa 2 {^ /asX^ fat&)] 



11) 



Note that this kinetic equation has a formal stationary Eq. JSJ by taking non-interacting two-particle Floquet 

"Fermi-Dirac" solution f a = [e^( Aa_ w + However, wave functions to be direct products of one-particle ones, 

this (3 has nothing to do with the real electronic temper- It has the obvious dependence on the positions: 
ature, which for this solution is infinite since the states 

with all s participate equally. Instead, f a will be as- Af ct3Q , 4 _ >QlCt2 ~ Al\ 3 \ i ^\ 1 \ 2 e~\ Sl+S2 ~ S3 ~ S4, \^ , (12) 

sumed to be independent of A, which corresponds to the 

actual distribution function f e to be smooth on the scale which comes from the overlap of the localized wave func- 

of single-particle energies e ~ T», while the dependence tions, the preexponential factor representing the matrix 

of f a on s determines f e through e = sui. element for \s\ + s 2 — s 3 — s 4 | ~ £. To estimate it one 

The matrix element M a3ai _+ aia2 originating from a can use the fact that the particle motion within one lo- 

two-particle interaction C/(r!,r 2 ) can be estimated from calization segment is diffusive. The starting point is 



\Mx 3 \ 4 ^\ 1 \ 2 



EE/ d ^ ■■■ d d r4xt 1 \r 1 ) X %\r 2 )}*U(r 1 ,r 2 )x^ (r 2 ) X ^ + * a -'»>(r 1 ) X 



Si ,S2 ,S 2 S3 ,54 ,S 



xx A s 1 3j (r3)x A : 4J (r 4 )t/(r3,r 4 ) ^(r 4 ) X^'^ (r 3 ) 



^r) E / U (^^2) (G R - G A )(r 2 , 4; r 4 , s' 4 ; A3) (G R - G A )(r 4l s 4 ; r 2 , s 2 ; A 2 ) 5 

xC/(r 3 , r 4 ) {G R - G A )(r 3 , s 3 ; n, Sl ; A x ) (G R - G A )(r uSl + s 2 - s' 2 ; r 3 , s 3 + s 4 - s 4 ; A 4 ), 

I 



the retarded and advanced Green's functions of the 
Schrodinger equation (JHJ) introduced in the usual way: 

xl' } (r) [x^ (*')}* 



G R ' A (r, S ;r',s'- 1 X) = J2- 



A - A' ±io 



(13) 



where the quasienergy A' labels the states within the seg- 
ment. It is technically simpler to consider the electron 
motion in the dot to be diffusive rather than ballistic 



(the final results are the same |lq). When averaging the 
products of Green's functions, one introduces the diffu- 
sion ladder dressed by the coupling V to the second or- 
der. The derivation, though cumbersome, is quite stan- 
dard 0| . The part of the screened Coulomb interaction, 
responsible for collisions, has the form (lif : 



t/(ri,r 2 ) 



L d 6 



S(ri - r 2 ) , 



(14) 
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and the final result is: the dot in the usual sense 

[2 5 4 ^ 1 S 4 



Ea?~«*7- (15) 



a 3 a 4 AiA 2 (27r) 2 £ A 2 C-E'Th Passing to the distribution function f e as mentioned 

above (integration over A2 , A3 , A4 in the collision integral 
Here Aj are the eigenvalues of the diffusion operator in gives a factor cj 2 /<5|), one arrives at the kinetic equation: 



dt 



St[/] = J ^^^lw(e 1 + e 2 -e3-e4)[(l-/ ei )(l-/ £2 )/e 3 /e 4 -/ ei /e 2 (l-/e3)(l-/e 4 )], (16) 



W(fi) 



-|C|/T, 



w(n) = w(-n), (17) 



the last property following from the isotropy of the spec- 
trum. If W(fl) were a (^-function, the collision inte- 
gral would be zero for a Fermi-Dirac distribution / e T = 
[e e / T + If the characteristic temperature T » T*, 

it is almost a (5-function, so is a good approximation. 
However, its temperature can change with time, so the 
solution should be sought in the form: 



dt 



dT dfl 



dt dT 



dt 



(18) 



where f e is small in T*/T. St[/] can be linearized with 
respect to /: St[/ e T + f e ] w St[/ £ T ] + St / e . In St one can 
approximate W(fi) oc 5(Q), then St(9/ e T /<9T) = 0, so 
/ e should be taken orthogonal to df^ /dT (as a function 
of e). The orthogonal component of the kinetic equation 
determines then the single-particle relaxation rate of the 
Floquet states due to collisions: 



Jirpl 



le(T) = 



2d 2 



2 rji'2 



TV T 



E Th 



(19) 



This coincides with the SIA expression |9J, T* having 
dropped out ^(| . The component of the kinetic equation 
along df T /dT determines the heating: 



dTdf? „ 15 -tv 2 dfl 
dt dT ~ 4(tt 2 - 6) dT 



W(fi) 



il 2 dfl 



(20) 



Estimating the integral, one arrives at Eq. (JJ. As a re- 
sult, the temperature is growing exponentially as given 
by Eq. (JIJ as long as 7 e (T) <C 1/t*, otherwise Fermi 
Golden Rule Q is no longer applicable, the Floquet 
states are destroyed, and the absorption becomes Ohmic. 

In conclusion, electron-electron collisions were studied 
in the regime of the strong dynamic localization in an 
ac driven quantum dot. It was shown that in this case 
the total electronic energy is effectively not conserved in a 
collision, which leads to heating of the electrons in the dot 
described by the kinetic equation l|16|) . As the frequency 



of collisions increases with temperature, this heating is 
exponential in time, as given by Eq. This picture 
persists until the strong dynamic localization is destroyed 
and the system returns in the Ohmic absorption regime. 

The author is greatly indebted to V. E. Kravtsov, 
D. L. Maslov, B. N. Narozhny, and O. M. Yevtushcnko 
for helpful discussions and a critical reading of the paper. 
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